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摘 要 : 以 玉龙 喀什 河流 域 为 研究 区 ,使 用 中 国 气象 同化 驱动 数据 集 (the China Meteorological Assimilation Driving 
Datasets for the SWAT Model, CMADS) 驱 动 SWAT(Soil and Water Assessment Tool) 水 文 模 型 进行 径流 模拟 ,评估 了 
CMADS 在 该 流域 的 精度 以 及 对 SWAT 模 型 的 适用 性 。 结 果 表 明 :(1) CMADS 的 降水 .最 高 气温 和 最 低 气 温 与 地 面 
水 文 测 站 获取 的 相应 参数 的 相关 系数 分 别 达到 了 0.650 .0.998 和 0.995 ,对 气温 的 模拟 结果 明显 优 于 降水 ;并 利用 协 


同 克 里 金 插值 分 析 了 CMADS 在 研究 区 的 降水 、 最 高 气温 和 最 低 气温 的 时 空 分 布 特征 ,表明 该 数据 集 能 够 很 好 地 反 


映 下 垫 面 地 形 特征 ,具有 较 高 的 质量 。(2) SWAT 模 型 在 模拟 月 尺度 径流 方面 , 除 极 值 模拟 较 差 以 外 , 率 定期 和 验证 
其 的 纳什 效率 系数 (Nash-Sutcliffe Efficiency Coefficient, NSE) 分 别 达到 了 0.845 和 0.836, 取 得 了 较为 满意 的 模拟 结 


果 。 本 文 证 明了 CMADS+SWAT 模 式 对 高 寒山 区 的 水 文 模拟 有 很 好 的 适用 性 ,提供 了 在 地 面 站 点 稀缺 地 区 构建 水 
文 模型 的 蔡 代 方案 ,对 寒 旱 区 的 水 文 模拟 起 到 了 一 定 的 促进 作用 。 
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玉龙 喀什 河 位 于 和 田 流 域 , 属 于 西北 内 陆 干 时 
区 ,水 资源 严重 制约 着 该 区 域 的 经 济 发 展 ,在 近 50 a 
来 和 田 流域 的 气候 特征 有 了 明显 的 变化 ,洪水 和 干 
旱 等 极端 气候 频繁 发 生 , 严 重 影响 了 人 类 活动 ,最 
需 精确 的 水 文 模型 来 模拟 预测 该 流域 的 地 表 径流 ， 
为 水 资源 和 洪水 风险 管理 提供 科学 指导 。 通 常 ， 
水 文 模型 采用 地 面 气象 站 收集 的 气象 数据 来 进行 
水 文 过程 和 径流 预测 的 研究 。 但 由 于 环境 特殊 、 
技术 条 件 落后 等 因素 ,在 西北 内 陆 干旱 区 ,尤其 是 
高 寒山 区 地 面 气象 测 站 稀 玻 ,甚至 无 地 面 气象 站 ， 
即使 有 测 站 的 地 方 也 存在 不 同 程度 的 缺 测 记录 ”。 
因此 ,地 面 气象 数据 的 不 完备 严重 制约 着 水 文 模 型 
的 模拟 精度 ,迫切 需要 更 高 分 辩 率 的 气象 数据 驱动 
水 文 模型 。 

目前 ,有 2 种 高 分 辩 率 气象 数据 :一 种 是 全 球 气 
候 模 式 或 区 域 气 候 模 式 的 降 尺 度数 据 , 男 一 种 是 大 
气 再 分 析 数 据 集 “。 相 比 于 前 一 种 ,大 气 再 分 析 数 
据 由 于 其 覆盖 范围 广 和 更 高 分 辩 率 的 优点 成 为 驱 
动 水 文 模型 的 首选 ,被 国内 外 学 者 广泛 使 用 于 各 类 
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水 文 模型 模拟 研究 "。 大 气 再 分 析 数 据 是 由 地 面 实 
测 数 据 和 气象 卫星 数据 再 结合 的 气候 模式 ,在 真实 
的 环流 场 强迫 下 ,通过 数据 同化 等 技术 制作 得 到 历 
史 气 象 数据 集 ”"。 国 际 上 使 用 较为 广泛 的 主要 有 美 
E] NCEP 气候 再 分 析 数 据 [Climate Forecast System 
Reanalysis (CFSR) , National Center for Atmospheric 
Research (NCAR) ] .美国 国家 航空 航天 局 再 分 析 产 
i (Modern-Era Retrospective Analysis for Research 
and Applications , MERRA ) .欧洲 中 期 天 气 预报 中 心 
(European Centre for Medium-Range Weather Fore- 
cast, ECMWF ) 的 EAR-Inter-Im ,EAR- 15, EAR-40 和 
中 国 大 气 同化 数据 集 CMADS 等 %*"。 许 多 学 者 对 
CMADS 在 中 国 不 同 地 区 的 适用 性 进行 了 验证 ,在 中 
国 西部 地 区 , 孟 现 勇 等 以 黑河 流域 为 例 , 使 用 
CMADS 、CFSR 和 TWS 再 分 析 数 据 驱动 SWAT 水 文 
模型 ,系统 地 分 析 了 3 类 再 分 析 数 据 对 黑河 流域 水 
文 过 程 的 影响 ,结果 表明 使 用 CMADS 再 分 析 数 据 
集 的 径流 模拟 效果 均 优 于 其 他 2 个 ; 张 利 敏 等 ”使 
用 CMADS 驱动 SWAT 模 型 研究 了 其 在 辽宁 浑 河流 
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域 的 适用 性 ,认为 CMADS 能 够 很 好 地 表现 下 热 面 
地 形 特征 ;Jun 等 和 刘 兆 晨 等 中 对 比分 析 了 
CMADS 和 CFSR 的 气候 要 素 与 气象 站 的 观测 值 ,并 
结合 SWAT 模 型 模拟 水 文 过 程 ,研究 表明 CMADS+ 
SWAT 模 式 的 模拟 效果 最 优 。 通 过 以 上 多 地 区 ,多 
方面 的 对 比 应 用 研究 ,基于 卫星 和 再 分 析 降 水 数据 
的 广泛 应 用 , 极 大 地 促进 了 地 面 实测 气象 站 点 稀少 
地 区 的 水 文 气象 研究 , 且 CMADS 再 分 析 数 据 集 比 
其 他 再 分 析 数 据 集 表现 出 更 高 模拟 潜力 。 

玉龙 喀什 河流 域 地 形 复杂 ,属于 典型 的 高 寒山 
区 ,流域 内 无 气象 站 点 分 布 ,不 能 满足 分 布 式 水 文 
模型 模拟 的 气象 驱动 要 求 ,必须 利用 气候 模式 或 再 
分 析 资 料 作为 流域 内 气象 替代 数据 。CMADS 作为 
适应 于 分 布 式 水 文 模型 的 气象 驱动 数据 ,还 尚未 在 
中 国 西北 寒 旱 区 玉龙 喀什 河流 域 得 到 充分 的 研究 ， 
需要 评估 其 在 研究 区 模拟 气温 和 降水 的 能 力 ,为 玉 
龙 喀什 河流 域 水 文 模型 的 构建 提供 备 选 数据 集 ,对 
研究 该 流域 水 资源 的 形成 和 分 布 规律 ,理解 高 寒山 
区 水 文 循环 过 程 具有 重要 科学 意义 ,对 促进 玉龙 喀 
什 河流 域 水 资源 合理 利用 和 防范 夏 洪 春 旱 等 自然 
灾害 提供 一 定 科学 指导 。 


1 数据 与 方法 


1.1 研究 区 概况 
玉龙 喀什 河 (以 下 简称 玉河 ) 为 和 田 河 源流 ( 东 
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支 ) 之 一 ,位 于 塔里木 盆地 南 缘 ,发 源 于 昆仑 山北 
坡 , 自 南 向 北 流 经 和 田 绿 洲 , 经 出 山口 同 古 孜 洛克 
水 文 站 (79.91609?°E,36.81055°N) 与 喀 拉 喀什 河 ( 西 
支 ) 并 行 , 在 下 游 于 阔 什 拉 什 汇流 入 和 田 河 ,最 终 注 
人 塔里木 河 ,全 长 504 km, WRK T E JEX 
什 河 上 游 河源 产 流 区 (77.25° ~81.75° E, 34.75? ~ 
36.25*N ,图 1) ,地 势 南 高 北 低 ,落差 极 大 , 集 水 面积 
为 14887 km ;研究 区 属于 内 陆 干旱 区 ,气候 干燥 、 降 
水 稀少 ,蒸发 强烈 中, 玉河 出 山口 同 古 孜 洛克 水 文 
站 日 均 气 温 12.7 % ,年 均 降水 94.2 mm ”。 玉 河 地 
表 径 流 补给 主要 依靠 冰川 积 雪 融 水 及 部 分 高 山 降 
水 ,平均 年 径流 量 为 22.6x10* m, FR ENEK, 
6 一 9 月 为 玉河 洪水 期 ,其 地 表 径 流量 占 到 全 年 的 
80% ,最 大 年 径流 超出 平均 径流 64.2% ,最 小 年 径流 
低 于 平均 径流 34.9%。 

1.2 数据 来 源 与 处 理 

12. 气象 数据 气象 数据 作为 分 布 式 水 文 模 型 的 
主要 驱动 数据 ,其 精确 性 与 完整 性 至 关 重 要 ”。 由 
于 研究 区 位 于 高 寒山 区 ,地 形 复杂 多 变 ,人 研究 区 内 
没有 相应 的 气象 观测 站 ,只 有 同 古 孜 洛克 水 文 站 ， 
但 该 站 位 于 流域 边缘 平原 处 ,无 法 替代 整个 流域 的 
气温 降水 情况 , 故 本 文采 用 中 国 大 气 同 化 驱动 数据 
集 CMADS V1.1 版 本 作为 SWAT 模 型 的 气象 驱动 数 
JE, CMADS 数据 集 数据 源 为 国家 级 自动 气象 站 
观测 数据 ,包括 测 站 的 经 纬度 和 海拔 高 度 以 及 降 
水 .气温 气压. 比 湿 、 辐 射 和 风速 等 气象 要 素数 据 ， 
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注 : 底 图 采用 自然 资源 部 标准 地 图 制作 , 审 图 号 为 GS(2019)3333 号 ,对 底 图 边界 无 修改 。 下 同 。 
图 1 研究 区 地 形 及 站 点 分 布 
Fig. 1 Topography and station distribution of the study area 
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融合 卫星 数据 产品 ,利用 多 重 网 格 三 维 变 分 ,循环 
般 套 、 重 采样 等 方法 ,并 以 NCEP/GFS 再 分 析 数 据 为 
背景 场 做 地 形 调 整 并 插值 到 0.25°x0.25° 格 点 上 。 
CMADS 数 据 集 已 按照 SWAT 模 型 气象 数据 输入 格 
式 进行 了 调整 ,不 需要 建立 每 个 站 点 的 天 气 发 生 融 
数据 库 , 只 需 建 立 研 究 区 的 气象 站 索引 表 。 研 究 区 
内 CMADS 站 点 分 布 见 图 1, 数 据 来 源 见 表 1。 

1.2.2 下 热 面 数据 ”SWAT 模型 所 需 下 热 面 数据 包 
FE DEM 数字 高 程 、 土 地 利用 数据 、 土 壤 数 据 ( 表 1)。 
DEM 数字 高 程 是 获取 研究 区 地 形 地 貌 .坡度 的 重要 
依据 ,用 以 生成 河 网 水 系 、 划 分 子 流域 及 确定 流域 
边界 2 ,经 过 SWAT 模 型 处 理 牛 成 的 玉河 水 系 以 及 
流域 边界 (图 1); 土 地 利用 数据 是 反映 研究 区 某 个 
时 期 的 土地 表面 要 素 的 空间 分 布 状态 、 地 表 特 征 以 及 
人 类 对 土地 的 开发 利用 的 数据 资料 中 ,对 于 土地 利 
用 数据 , 需 进行 土地 利用 重 分 类 预 处 理 , 之 后 与 
SWAT 模 型 内 部 数据 库 的 索引 表 进 行 链接 。 根 据 中 


国 科学 院 资源 环境 科学 数据 中 心 分 类 指标 和 研究 
区 实际 情况 ,将 玉河 土地 利用 类 型 重 分 类 为 6 类 (图 
2) ,分 别 是 耕地 、 林 地 .草地 水域 .永久 性 冰川 雪 
地 、 裸 地 ,其 中 裸 地 的 占 比 最 大 ,为 49.23% ;土壤 数 
据 用 以 反映 研究 区 不 同 土壤 的 空间 分 布 以 及 各 种 
土壤 物理 属性 ;通过 对 世界 土壤 数据 库 的 研究 区 范 
围 裁剪 以 及 重 分 类 后 ,得 到 研究 区 的 土壤 类 型 空间 
分 布 ( 图 3), 其 中 永 冻 薄 层 土 占 比 最 大 ,为 36.04%; 
对 于 土壤 物理 属性 ,部 分 土壤 基本 数据 可 以 从 
HWSD 土壤 数据 库 获 得 ,如 土壤 各 粒 径 成 分 的 含量 、 
碳酸 盐 含 量 ` 有 机 碳 含量 等 ,大 部 分 数据 无 法 直接 
测 得 , 需 通 过 SPAW 软件 的 Soil Water Characteristic 
模块 计算 参数 ,如 土壤 湿 密 度 、 土 壤 有 效 含水 量 、 饱 
和 水 力 传导 系数 等 。 

1.3 研究 方法 

1.3.1 SWAT 模 型 设置 SWAT 模 型 是 由 美国 农业 
部 农业 研究 中 心 开 发 的 基于 流域 尺度 的 一 个 长 时 


表 1 数据 类 型 及 来 源 


Tab.1 Data types and sources 


数据 名 称 数据 描述 数据 来 源 
气象 数据 zx [h] 4] 9E 0.25? ,时 间 为 2008 一 2016 年 逐日 国家 青藏 高 原 科 学 数据 中 心 ,数据 版 本 CMADS V1.1 
水 文 数据 1969 一 2016 年 逐 月 径流 和 田地 区 水 利 局 ,玉龙 喀什 河 出 山口 同 古 孜 洛克 水 文 站 
DEM 数字 高 空间 分 辩 率 30 m, 地 理 坐 标 为 WGS_1984, 投影 坐 标 为 ”地 理 空间 数据 云 ASTER GDEM 30 m 分 辨 率 数 字 高 程 数 据 
程 数据 WGS. 1984. World. Mercator 
土地 利用 数据 “空间 分 辨 率 30 m, 时 间 为 2010 年 ,投影 坐标 为 WGS。 中国 科学 院 资源 环境 科学 数据 中 心 

1984 World. Mercator 

土壤 数据 空间 分 辩 率 1 km, 时 间 为 2008 4E ,投影 坐标 为 WGS_ 世界 粮农 组 织 提供 世界 土壤 数据 库 (HWSD) 中 的 中 国土 


1984 World Mercator 
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Fig. 2 Land use types 
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段 分 布 式 流域 水 文 模型 ,可 以 直接 模拟 水 流 、 泥 沙 
等 物理 过 程 ”, 可 以 定量 化 的 输入 参数 (如 气候 、 植 
被 等 变化 ) 以 研究 对 径流 等 的 影响 下 。 根 据 人 研究 
区 实际 情况 ,设置 流域 面积 划分 阀 值 为 25000 hm’, 
将 流域 划分 为 53 个子 流域 , 子 流域 平均 面积 > 
280.88 km ;将 流域 划分 为 1863 个 水 文 响应 单元 
HRU (Hydrologic Response Unit) 。 在 加 载 2008— 
2016 ^F. CMADS 气象 驱动 数据 之 后 ,建立 SWAT 月 
径流 模型 ,设置 2008 一 2009 年 为 预 热 期 ,2010 一 
2012 年 为 模型 参数 率 定 期 ,2013 一 2016 年 为 模型 验 
证 期 。 

模型 参数 的 率 定 采 用 SWAT-CUP(SWAT Cali- 
bration and Uncertainty Program ) 软 件 , 此 软件 专门 用 
于 SWAT 模 型 的 参数 敏感 性 分 析 、 模 型 的 率 定 以 及 
不 确定 性 分 析 “”1。 通 过 比较 SWAT-CUP 所 提供 的 
5 种 算法 ,SUFI-2 算 法 结合 了 不 确定 性 分 析 与 优化 ， 
并 使 用 拉丁 超 立 方 抽样 方法 处 理 参数 ,具有 运行 简单 、 
计算 效率 高 以 及 可 灵活 选择 目标 函数 等 优点 汪 ”。 
故 本 文选 用 SWAT-CUP 提 供 的 SUFI-2 算 法 ,并 选择 
纳什 效率 系数 (Nash-Sutcliffe Efficiency Coefficient， 
NSE) 为 目标 函数 来 率 定 模 型 参数 ,利用 i-Stat 值 和 


P-Value 值 进行 模型 参数 敏感 性 检验 。 

1.3.2 CMADS 数据 集 评价 方法 为 评价 CMADS 数 
据 集 在 研究 区 的 模拟 精度 ,选取 距离 同 古 孜 洛 
克 水 文 站 最 近 的 CMADS 的 148-81(80.03125?°E， 
36.78125°N) 站 点 进行 分 析 。 因 同 古 孜 洛克 水 文 站 
仅 有 逐日 降水 以 及 日 平均 气温 , 故 只 评价 CMADS 
数据 集 与 观测 站 点 的 降水 与 平均 气温 。 对 于 平均 
气温 主要 采用 相关 系数 、 均 方 根 误差 和 相对 误差 3 
项 指标 ;而 对 于 降水 为 了 全 面 表征 CMADS 数据 集 
在 该 地 区 模拟 的 效果 ,增加 了 探测 率 、 误 报 率 和 临 
界 成 功 指数 3 项 指标 ,具体 公式 如 表 2。 

1.3.3 径流 模拟 效果 评价 方法 ”对 SWAT 模 型 的 径 
流 模 拟 精 度 采 用 纳什 效率 系数 .模型 决定 系数 和 相 
对 误差 3 类 指标 进行 评价 ,具体 公式 如 表 2。 


2 结果 与 分 析 


24. CMADS 精度 及 时 空 分 布 特征 评价 

2.1.1 CMADS 精度 及 评价 ”降水 数据 的 精度 很 大 
程度 上 决定 了 模型 的 径流 模拟 效果 ,为 直观 分 析 
CMADS 降水 数据 与 研究 区 同 古 孜 洛克 水 文 站 实测 
降水 的 拟 合 程度 , 作 日 尺度 降水 过 程 线 (图 4)。 由 


表 2 精度 评价 公式 


Tab.2 Accuracy evaluation formula 


名 称 公式 TER 最 优 值 
相关 系数 (Correlation Coefficient, CC ) Y. O- V...) - V..)] [71,1] 1 
obs obs simV 7 — " sim 
CC= 一 
地 a X S q^ 
ÈO- V.P 31V. - V..GT 
1 Ni=1 
均 方 根 误 差 (Root Mean Square Error, RMSE) RMSE- DM D" [0,99] 0 
A Iis 
相对 误差 (Relative Error, RE ) bU -VOl [799,9] 0 
sim obs. 
RE=5 x 100% 
ÈVO 
探测 率 (Probability of Detection, POD ) POD- ..H [0.1] 1 
H+M 
误 报 率 (False Alarm Ratio, FAR ) FAR= F [0.1] 0 
H«F 
临界 成 功 指数 (Critical Suecess Index, CSI) CSI= H [0,1] 1 
H+M+F 
纳什 效率 系数 (Nash-Sutcliffe Efficiency Coefficient, NSE) Y à- V.P [0.1] 1 
obs\%7 sim 
NSE =1- -+= 


izl 


— me m t . 2 
模型 决定 系数 (Coefficient of Determination, R°) or. -7 VO- UE [0.1] 1 
Rc 


PO FS Va- V. 


注 :V.( 让 表示 实测 变量 [具体 可 表示 为 降水 (mm) ,气温 (C), 径 流量 (mm )]; Vi,(i) 表 示 模 拟 变 量 ;n 表示 模拟 数 ; 太 表 示 模 拟 与 地 面 观测 均 有 降 
水 的 天 数 (d); 必 表示 模拟 无 降水 而 地 面 有 降水 的 天 数 (d) ;Ff 表示 模拟 有 降水 而 地 面 无 降水 的 天 数 (d)。 
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图 4 同 古 孜 洛克 水 文 站 与 CMADS 148-81 
格 点 日 降水 过 程 对 比 
Fig.4 Comparison of precipitation process of Tongguziluoke 
hydrological station between and grid 148-81 of CMADS 


图 4 可 以 看 出 ,CMADS 降 水 所 表现 出 的 峰值 以 及 峰 
值 出 现 的 时 间 与 观测 降水 在 大 部 分 年 份 都 有 较 好 
的 对 应 ,但 在 个 别 年 份 存在 峰值 模拟 效果 较 差 ,如 
在 2010—2011 年 出 现 明显 的 错 峰 现象 ,在 2013 一 
2014 年 CMADS 的 降水 峰值 明显 低 于 实测 降水 ,说 
明 CMADS 对 连续 的 极端 降水 和 变化 剧烈 的 降水 事 
件 的 模拟 性 能 较 差 。 

为 定量 评估 CMADS 降水 数据 的 拟 合 效果 , 作 
散 点 图 进行 线性 拟 合 ( 图 5)。CMADS 降 水 与 实测 
降水 的 相关 系数 CC 为 0.65, 两 者 之 间 具 有 较 强 的 线 
性 相关 性 ;相对 误差 RE 为 24.67% 大 于 0,CMADS 降 
水 与 实测 降水 数据 年 均 降 水 量 分 别 为 97.77 mm 与 
91.03 m, CMADS 对 实测 降水 有 所 高 佑 ,但 相差 仅 
6.74 mm; 均 方 根 误差 RMSE 为 1.14, 表 明 两 者 之 间 
误差 较 小 ;探测 率 POD RIK FAR 和 临界 成 功 指 


35 
CC-0.65 o 实测 降水 与 CMADS 的 散 点 

30 L Hnc ---- 最 高 气温 的 线性 拟 合 
Pa -1.14 一 jy=x 参 考 线 。 
T 25| POD-087 ^ ^. 
S FAR-0.25 
E 20 上 CSI-0.68 
SA 6 
a 15 2 
5 10 p 

[P 
5 上 2e i94 o 
d go $ o 
0 o 9 o 


OQ | L L 1 L 


实测 /mm.d) 


图 5 同 古 孜 洛 克 水 文 站 与 CMADS 的 
148-81 格 点 降水 散 点 图 
Fig. 5 Precipitation scatter diagram of Tongguziluoke 
hydrological station between and grid 148-81 of CMADS 


数 CSI 分 别 为 0.87.0.25 和 0.68 ,说 明 CMADS 对 于 实 
际 降水 事件 发 生 的 概率 拟 合 准确 性 较 高 。 总 体 而 
言 ,CMADS 降水 数据 不 论 在 降水 量 还 是 降水 事件 再 
现 的 模拟 上 ,都 有 着 较 好 的 表现 ,CMADS 降水 数据 
在 研究 区 有 较 高 的 精度 和 可 靠 性 ,加 上 其 覆盖 全 
面 ,分 辩 率 高 的 特点 ,可 在 研究 区 代替 实测 降水 数 
据 驱 动 SWAT 水 文 模型 进行 径流 的 模拟 。 

气温 作为 SWAT 模 型 必 不 可 少 的 气象 数据 输入 
项 ,决定 着 SWAT 模 型 降水 形态 的 转变 以 及 融雪 速 
度 和 融雪 水 量 。 人 研究 区 属于 内 陆 高 寒山 区 ,融雪 是 
径流 的 重要 组 成 部 分 , 故 对 融雪 模拟 的 效果 也 将 对 
最 终 径 流 的 模拟 效果 产生 很 大 的 影响 。 相 对 于 降 
水 来 说 ,CMADS 数据 集 对 气温 的 拟 合 效 果 要 优 于 降 
水 。 从 气温 拟 合 散 点 图 (图 6 和 图 7) 可 知 ,对 于 最 高 
气温 和 最 低 气 温 的 拟 合 ,CMADS 气 温 与 实测 最 高 气 
温 与 最 低 气 温 之 间 的 相关 系数 CC 分 别 达 到 了 0.998 
与 0.995,CMADS 和 气温 与 实测 气温 之 间 具 有 相当 高 
的 线性 相关 性 ;相对 误差 RE 分别 为 11.22% 与 
16.52% 均 大 于 零 ,CMADS 对 实际 气温 是 低估 的 ; 均 
方 根 误差 RMSE 分 别 为 1.787 和 1.75,CMADS 对 实 
际 气温 的 模拟 误差 较 小 。 通 过 以 上 指标 分 析 ， 
CMADS 对 于 气温 的 模拟 具有 相当 高 的 精度 和 可 靠 
性 ,可 用 于 研究 区 作为 SWAT 模 型 的 气象 驱动 数据 
来 模拟 径流 。 
2.1.2 CMADS 降水 气温 时 空 分 布 特征 评价 SWAT 
模型 作为 分 布 式 水 文 模型 ,不 仅 对 输入 数据 精度 有 


CC=0.998 
30 F RE=11.22% 
RMSE=1.787 dh 
| o 实测 最 高 气温 与 CMADS M 


A 最 高 气温 的 线性 拟 合 
c 10 
9 o 
-10 
-20 
-15 -10 0 10 20 30 35 
实测 /*C 


图 6 同 古 孜 洛克 水 文 站 与 CMADS 的 148-81 
格 点 最 高 气温 散 点 图 
Fig.6 Maximum temperature scatter diagram of 
Tongguziluoke hydrological station between and 
grid 148-81 of CMADS 
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图 7 同 古 孜 洛克 水 文 站 与 CMADS 的 148-81 
格 点 最 低 气 温 散 点 图 
Fig.7 Minimum temperature scatter diagram of 
Tongguziluoke hydrological station between 
and grid 148-81 of CMADS 


较 高 的 要 求 ,而 且 对 气温 与 降水 的 空间 分 布 特征 也 
有 和 较 高 的 要 求 。 从 地 形 图 (图 1) 可 以 看 出 ,研究 区 
地 形 复杂 ,海拔 变化 大 ,36°15'N 为 流域 的 高 程 划 分 
线 , 其 下 至 流域 出 水 口 , 是 山 前 冲积 扇 平 原 区 ,地 势 
平坦 ;其 上 至 流域 边界 ,海拔 均 在 3000 m 以 上 ,是 冰 
川 积 雪 的 主要 覆盖 地 , 沟 殉 交错 。 为 充分 考虑 地 形 
对 降水 和 气温 的 影响 ,采用 ArcGIS 中 地 统计 处 理 模 
块 的 协同 克 里 金 法 对 CMADS 的 气温 与 降水 进行 空 
间 插 值 ,协同 克 里 金 插值 是 基于 高 程 利 用 区 域 化 变 
量 的 原始 数据 和 变异 函数 的 结构 特点 ,对 未 采样 点 
的 区 域 化 变量 的 取 值 进行 线性 无 偏 、 最 优 估计 。 

对 于 降水 ,从 时 间 分 布 上 来 看 (图 8) ,年 内 降水 
主要 集中 于 夏季 ,夏季 流域 降水 量 占 全 年 总 降水 量 
的 38.0396 , 春季 降水 最 少 , 占 全 年 总 降水 量 的 
18.52% ;从 空间 分 布 来 看 ,降水 并 未 表现 出 明显 的 
随 海拔 升 高 而 降水 量 增加 的 趋势 ,四 季 的 结果 均 表 
现 出 强烈 的 地 形 特征 的 降水 集中 区 ,图 8 中 所 表现 
的 降水 集中 区 刚好 位 于 流域 内 一 处 开口 向 东 的 马 
蹄 形 盆地 ,气流 从 东 面 开口 进入 ,使 气流 辐 合 上 升 ， 
形成 地 形 雨 。 此 外 ,四 季 均 存在 一 条 20~30 mm 的 
降水 分 割 带 ,该 降水 带 主要 沿 昆仑 山脉 走势 分 布 。 
以 北 由 于 昆仑 山 阻 挡 了 南 来 的 暖 空气 ,造成 高 山区 
降水 较 少 ;同时 ,该 降水 带 的 移动 也 表征 了 四 季 的 
变换 。 总 体 而 言 ,CMADS 降水 数据 的 空间 分 布 能 够 
较 好 地 反映 研究 区 的 地 形 特征 。 

对 于 气温 ,从 时 间 分 布 来 看 (图 8) ,日 最 高 气温 
均值 出 现在 秋季 ,达到 28 % ,日 最 低 气 温 均值 出 现 


在 春季 和 冬季 ,达到 -28 % ;以 日 最 高 气温 均值 为 
例 , 自 春 季 开 始 ,4~8 Y 气 温带 逐渐 向 南 转移 ,秋季 
到 达 流 域 中 部 ,之 后 再 次 转移 到 流域 最 北端 ,此 为 
年 内 的 气温 循环 。 从 空间 分 布 来 看 ,整体 上 呈现 出 
随 海拔 升 高 气温 降低 的 趋势 ,气温 直 减 率 约 为 
-1.56 .km ', 其 符合 高 寒山 区 的 气温 直 减 率 ”。 
另外 ,四 季 气 温 空 间 分 布 均 存在 一 个 气温 低 值 区 ， 
位 置 恰好 位 于 降水 集中 区 ,并 在 此 区 域 形 成 环 状 气 
温带 , 易 形成 气流 辐 合 ,从 男 一 方面 说 明了 该 区 域 
为 降水 集中 区 的 原因 。 通 过 以 上 对 降水 和 气温 的 
时 空 分 布 特征 分 析 , 说 明 CMADS 很 好 地 捕捉 到 了 
地 形 的 变化 。 
2.2 模型 参数 率 定 结果 

根据 研究 区 实际 情况 ,拟定 17 个 参数 进行 率 
定 , 因 人 研究 区 位 于 高 寒山 区 ,径流 补给 主要 为 高 山 
降水 .地 下 水 和 冰川 融雪 , 故 参 数 主要 根据 径流 补 
给 来 源 选择 ( 表 3),t-Stat 值 表示 敏感 性 的 程度 ,绝对 
值 越 大越 敏 感 ;P-Value 值 决定 敏感 性 的 显著 性 , 值 
越 接近 0, 越 显著 。SUFI-2 为 SWAT 模 型 最 常用 的 
参数 识别 算法 ,该 算法 识别 得 到 的 参数 最 优 值 可 能 
落 在 参数 先 验 分 布 范 围 之 外 , 当 模 型 率 定 出 现 不 合理 
参数 取 值 时 ,应 对 这 些 参数 取 值 进行 分 析 和 调整 3。 

在 敏感 性 排序 前 5 的 参数 中 ,与 融雪 模块 和 地 
下 水 模块 相关 的 参数 各 占 2 个 ,其 中 最 大 融雪 因子 
(SMFMX ,北半球 为 6 月 21 号 的 融雪 因子 ) 最 为 敏 
感 ,最 优 值 为 5.1345 mm - C^ - 4^ ; hz eil SMTMP 
为 0.5 %C, 当 子 流域 内 的 积 雪 温 度 超过 此 阀 值 将 发 
生 融 雪 , 以 上 参数 值 与 同属 于 寒 旱 区 的 开 都 河 月 尺 
度 径流 模拟 参数 取 值 相近 六 。 对 于 地 下 水 模块 , 基 
流 a 因 子 ALPHA_BF 表现 最 为 敏感 , 它 是 地 下 水 径 
流 对 补给 量变 化 响应 的 直接 指标 , 率 定 最 优 值 为 
0.1100; 其 次 为 地 下 水 延迟 时 间 GW_DELAY , 它 表 
征 了 地 下 水 流出 土壤 剖面 补给 地 表 水 的 时 间 延 迟 ， 
根据 尤 扬 等 “在 不 同 气候 情景 下 和 田 河 上 游 径流 
变化 中 ,建立 的 SWAT 月 尺度 模型 确定 参数 ,最 终 率 
定 为 329 d;TLAPS 最 终 率 定 结果 为 1.532 "C km! , 
与 气温 的 垂直 分 布 特征 所 表现 的 结果 一 致 。 通 过 
以 上 分 析 ,SWAT-CUP 率 定 的 各 项 参数 和 敏感 性 排 
序 ,与 刘 晓 稍 55 基 于 SWAT 模 型 的 和 田 河 上 游 气候 
和 土地 利用 变化 的 水 文 效 应 模拟 研究 中 ,所 构建 的 
月 尺度 模型 的 参数 结果 较为 相近 , 均 处 于 合理 的 物 
理 范 围 内 。 
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图 8 CMADS 降水 气温 时 空 分 布 
Fig. 8 Temporal and spatial distribution of CMADS precipitation and temperature 


2.3 径流 模拟 结果 及 适用 性 评价 

由 表 4 可 知 ,模型 在 率 定 期 和 验证 期 的 纳什 系 
数 分 别 达到 0.845 和 0.836 ,模型 决定 系数 分 别 达 到 
0.8562 和 0.8153 , 验证 期 此 2 项 指标 均 有 所 减 小 。 
率 定 期 和 验证 期 实测 月 径流 量 均值 较 模拟 均值 分 
别 高 出 7.464 m°- s 110.975 nis" ,相对 误差 较 率 
定期 增 大 为 -12.47% ,模拟 的 均值 偏 低 和 相对 误差 
为 负 , 均 说 明 径流 的 模拟 值 低估 了 实测 值 。 但 总 体 
上 来 看 , 率 定期 和 验证 期 的 各 项 指标 较为 可 观 , 表 
明 SWAT 模 型 在 该 流域 的 径流 模拟 中 具有 较 好 的 适 
应 性 。 

模拟 期 内 模拟 径流 与 实测 径流 曲线 趋势 基本 


一 致 (图 9) ,在 每 年 5 月 径流 逐渐 增 大 ,于 8 月 达到 
峰值 ,在 9 月 又 降 , 从 10 月 开始 到 次 年 4 月 期 间 均 为 
基 流 , 基 流 天 数 约 210 d。 局 部 来 看 ,存在 部 分 误差 
较 大 的 时 间 点 ,其 中 最 大 负 误 差 出 现在 2012 年 8 
月 ,模拟 值 低 于 实测 值 246.8 ma .s ,最 大 正 误 差 出 
现在 2013 年 9 月 ,模拟 值 高 于 实测 值 244.3 m^ s, 
最 小 误差 为 2010 年 9 月 ,误差 仅 为 0.3 ms, KAR 
分 模拟 值 处 于 实测 值 的 下 方 (图 9) ,与 评价 指标 中 
相对 误差 所 得 到 的 结果 一 致 , 旦 低估 现象 主要 出 现 
在 基 流 期 和 洪水 期 。 基 流 期 一 般 从 冬季 到 次 年 春 
季 末 ,雨季 已 过 , 且 气 温 已 经 降 至 0% 以 下 ,冰川 和 
积 雪 融 水 又 减 ,径流 主要 由 地 下 水 补给 。 但 研究 区 


202210.00025v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


TIRE SE : CMADS 在 玉龙 喀什 河 径流 模拟 中 的 适用 性 研究 


表 3 参数 率 定 结果 


Tab.3 Parameter calibration results 


序号 调 参 方式 ”参数 名 称 调 参 区 间 最 优 值 t-Stat P-Value 
1 V. SMFMX 最 大 融雪 因子 /mm * Cd) 1.5,6.9] 5.1345 54.7774 0.0000 
2 R ALPHA BF ” 基 流 a 因子 /d 0,1] 0.1100 -25.2418 0.0000 
3 V. SMTMP Tk c Acc —5,10] 0.5000 8.4776 0.0000 
4 V GW DELAY ”地 下 水 延迟 时 间 /d 0,500] 329.000 4.8323 0.0000 
5 A. TLAPS 气温 直 减 率 /(% km!) -1.5,0] -1.5320 2.4267 0.0154 
6 R OV N 坡 面 漫 流 的 曼 宁 系数 -0.2,0.2] —0.0090 2.1911 0.0286 
7 R CN2 初始 SCS 径 流 曲线 数 -0.2,0.2] -0.1000 -1.5313 0.1259 
8 R. SOL BD 土壤 湿 容 重 /(g.em-”) -0.2,0.2] -0.0278 -1.4564 0.1455 
9 R SOLK 饱和 渗透 系数 /(mmh') -0.5,0.5] 0.3000 -1.2161 0.2241 
10 A. PLAPS 降水 直 减 率 /(mm:km) 0,2] 1.6500 -1.1239 0.2612 
11 V. RCHRG DP. ”深层 含水 层 的 渗透 系数 0,1] 0.2480 0.7908 0.4292 
12 V. SFTMP KEC -5,5] -0.0216 0.6371 0.5242 
13 R_HRU_SLP 平均 比 降 -0.2,0.2] 0.0587 -0.5362 0.5919 
14 R SLSUBBSN ”平均 坡 长 /m -0.1.0.1] 0.0097 0.4710 0.6377 
15 V. GWQMN 发 生 回 归 流 所 需 的 浅 层 含水 量 的 水 位 阀 值 /mm 0,5000] 4750.0 0.4302 0.6671 
16 V. SMFMN 最 小 融雪 因子 /(mm* d) 1.5,6.9] 4.4167 -0.2597 0.7951 
17 V. ESCO 土壤 蒸发 补偿 因子 0,1] 0.4100 0.1570 0.8753 


注 :V 表 示 现 有 的 参数 值 将 被 给 定 的 值 取代 ;A 表 示 将 多 


400 


300 


月 均 流 量 /(m3.s 1!) 


200 


100 


台 定 的 值 加 上 现 有 的 值 ;R 表 示 将 现 有 的 参数 值 乘 以 (1+ 给 定 的 值 )。 


2010 2011 


表 4 径流 模拟 评价 结果 


Tab.4 Runoff simulation evaluation results 


1 
2014 


I 
2013 


2012 
年 份 


图 9 月 径流 模拟 结果 


Fig.9 Runoffsimulation results 


平均 值 /ms) —m Mii a un 
时 期 E 纳什 系数 PENE 相对 误差 

实测 。 ”模拟 系数 偏 低 。 
率 定期 93.678 86.214 0.845 0.8562 -10.93% 
验证 期 89.446 78.471 0.836 0.8153 —12.4796 


2015 


2016 2017 


se 对 于 地 下 水 的 监测 和 研究 较 少 , 缺 
资料 ,无 法 对 参数 进 和 


了 准确 的 率 定 ,导致 模型 中 
吉 果 不 够 理想 ,造成 模拟 期 基 流 


模拟 结果 表明 ,校准 期 和 验证 期 内 ,玉河 径流 
量 模拟 值 和 实测 值 之 间 的 各 项 评价 指标 的 结果 与 
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型 在 西北 寒 旱 区 建立 的 月 尺度 径流 模拟 结果 相近 ， 
均 在 允许 范围 之 内 。 因 此 ,SWAT 模 型 总 体 上 适用 
于 该 区 域 高 寒山 区 的 水 文 过 程 模 拟 。 


3 讨论 


本 文 的 研究 结果 表明 :气象 数据 作为 水 文 模型 
的 输入 数据 ,其 不 确定 性 对 径流 模拟 结果 有 着 重要 
的 影响 号 。 由 于 研究 区 的 资料 限制 ,只 能 使 用 仅 有 
的 出 山 日 控制 站 同 古 孜 洛克 水 文 站 来 评价 CMADS 
在 此 位 置 的 精度 ,具有 较 大 的 不 确定 性 。 对 于 降 
水 ,其 不 确定 性 主要 在 于 单 站 点 的 降水 控制 面积 
限 ,无 法 代表 整个 流域 的 实际 降水 情况 ,包括 降水 
强度 变化 .降水 位 置 以 及 降水 事件 的 发 生 3 ,降水 
站 点 密度 对 水 文 建 模具 有 一 定 影响 ;对 于 气温 ， 
研究 区 海拔 变化 幅度 大 , 且 存 在 冰川 积 雪 ,实际 气 
温 的 变化 无 法 得 到 测量 ,所 以 无 法 准确 评估 CMADS 
降水 和 气温 时 空 变异 性 。 在 评价 CMADS 的 时 空 分 
布 特征 时 ,使 用 协同 克 里 金 插值 对 未 采样 点 的 区 域 
进行 取 值 ,得 到 的 结果 虽然 是 一 个 最 优 估计 值 ,但 
并 不 能 反映 未 采样 点 的 实际 值 ,其 空间 分 布 特征 也 
存在 一 定 误差 ,可 进一步 引入 海拔 、 坡 向 、 坡 度 等 因 
子 ,提高 插值 精度 ””。 整 体 而 言 ,CMADS 在 精度 以 
及 时 空 特征 分 布 方面 均 有 着 良好 的 表现 ,用 以 驱动 
SWAT 模 型 模拟 研究 区 的 径流 ,其 结果 有 着 较 高 的 
可 靠 性 。 

本 文 使 用 2008 一 2012 的 实测 径流 率 定 SWAT 
模型 ,并 以 率 定 好 的 模型 验证 2013—2017 的 实测 径 
流 ,结果 表明 验证 期 的 指标 较 率 定期 的 差 ,模拟 效 
果 较 差 不 仅 是 因为 产 流 模式 与 模型 预 设 不 匹配 、 验 
证 期 水 文 环境 变化 等 ,还 与 率 定 径流 序列 过 短 导 致 
率 定期 模型 水 文 过 程 表 征 不 准确 有 关 ”。SWAT 水 
文 模型 作为 分 布 式 水 文 模型 ,尽管 在 结构 上 考虑 了 
融雪 和 冻 土 对 水 文 循环 的 影响 吕 ,拥有 较 多 的 物理 
参数 ,有 着 相对 完善 的 水 文 循环 过 程 , 但 是 由 于 缺 
乏 冰 川 模块 ,在 夏季 汛期 缺少 冰川 融 水 资料 , 故 造 
成 模拟 过 程 中 模拟 极 值 达 不 到 实测 极 值 。 孙 铭 悦 
等 加 使 用 格 点 数据 驱动 HBV 水 文 模型 在 西北 寒 旱 
区 呼 图 壁 河流 域 的 日 尺度 径流 模拟 中 ,同样 存在 
模拟 极 值 低 于 实测 径流 极 值 的 现象 。 祖 拜 代 : 木 
依 布 拉 等 “使 用 单 站 点 构建 了 乌鲁木齐 河上 游 的 


月 尺度 SWAT 模 型 ,验证 期 决定 系数 仅 为 0.75, 且 相 
对 误差 高 达 22.20% , 相 比 之 下 ,本 文 以 CMADS 格 点 
数据 驱动 SWAT 模 型 在 寒 旱 区 的 径流 模拟 效果 要 更 
好 , 故 CMADS+SWAT 模 型 在 寒 旱 区 玉河 流域 有 着 
较 好 的 适用 性 。 

4 结 论 

本 文 基于 研究 区 DEM 、2010 年 土地 利用 和 土壤 
数据 ,构建 了 玉龙 喀什 河 产 流 区 的 SWAT 模 型 ,验证 
了 CMADS 再 分 析 数 据 集 的 精度 和 空间 分 布 特征 ， 
模拟 了 研究 区 的 月 尺度 径流 ,得 出 如 下 主要 结论 : 

(1) 总 体 上 CMADS 再 分 析 数 据 的 精度 和 可 靠 
性 较 好 ,可 替代 无 站 点 流域 作为 气象 数据 驱动 水 文 
模型 。 气 温 的 模拟 效果 优 于 降水 ,降水 大 部 分 模拟 
值 与 实测 值 较为 接近 ,CMADS 再 分 析 数 据 的 时 空 分 
布 特征 基本 符合 实际 ,能 够 准确 地 捕捉 到 下 热 面 的 
地 形 特征 变化 ,但 CMADS 对 于 极端 的 气候 变化 的 
模拟 能 力 有 待 提高 。 

(2) 以 CMADS+SWAT 模 式 模拟 玉河 流域 月 尺 
度 径流 ,整体 上 模拟 径流 序列 与 实测 径流 变化 趋势 
一 致 ,在 研究 区 具有 较 好 的 适用 性 ,但 在 峰值 部 分 
出 现 明显 的 低估 现象 ,主要 原因 是 缺少 夏 汛 期 的 冰 
川 融 水 资料 ;另外 ,CMADS 数 据 集 本 身 的 时 间 跨 度 
不 够 长 ,导致 验证 期 的 模拟 效果 稍 差 ,模型 的 参数 
率 定 存在 一 定 的 偏差 。 

总 体 而 言 ,CMADS 再 分 析 数 据 集 通 过 研究 区 的 
验证 ,对 于 无 站 点 流域 的 水 文 模拟 效果 较 好 ,可 为 
无 降水 ,人 少 降水 站 点 或 缺失 降水 资料 的 区 域 提供 较 
为 可 靠 的 数据 来 源 , 为 扩大 水 文 模拟 的 时 间 和 空间 
尺度 提供 了 可 能 。 但 缺乏 与 其 他 再 分 析 数 据 的 比 
较 , 在 后 续 的 研究 中 将 耦合 其 他 再 分 析 数 据 进 行 比 
较 分 析 ,对 再 分 析 数据 进行 集合 平均 ,结合 各 类 再 
分 析 数 据 的 优势 ,减少 误差 ;并 收集 研究 区 冰川 资 
料 , 加 入 冰川 模块 ,完善 水 文 模型 结构 ,使 模拟 结果 
精度 更 高 。 
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Applicability of CMADS in runoff simulation of Yulong Kashi River 


LUO Chengyan, CHEN Fulong, HE Chaofei', 
LONG Aihua", QIAO Changlu' 
(1. College of Water & Architectural Engineering, Shihezi University, Shihezi 832000, Xinjiang, China; 2. State Key 
Laboratory of Simulation and Regulation of Water Cycle in River Basin, China Institute of Water Resources and 


Hydropower Research, Beijing 100038, China) 


Abstract: In this paper, the Yulong Kashi River Basin was selected as the research area, and the runoff simulation 
of SWAT hydrological model driven by The China Meteorological Assimilation Driving Datasets (CMADS) for 
the SWAT model was carried out to assess the accuracy of CMADS in the watershed and its applicability to the 
SWAT model. The results showed that (1) the correlation coefficients of precipitation, maximum temperature, and 
minimum temperature of CMADS with the measured hydrological stations were 0.650, 0.998, and 0.995, 
respectively, and their simulation of temperature was significantly better than that of precipitation. The spatial- 
temporal distribution characteristics of CMADS to precipitation and maximum and minimum temperatures were 
studied by using collaborative Kriging interpolation, which shows that the data set can well capture the 
topographic characteristics of the underlying surface in the study area, and its meteorological data has high 
quality. (2) With regard to runoff simulation of the SWAT model, except for poor simulation, the Nash-Sutcliffe 
efficiency coefficient of the periodic rate and validation period reached 0.845 and 0.836, respectively, and 
satisfactory simulation results were obtained. This paper showed that the CMADS + SWAT model had good 
applicability for hydrological simulation in the study area and provided an alternative to construct a hydrological 
model in areas with scarce ground stations, which plays a certain role in promoting hydrological simulation in 
cold and dry areas. 

Keywords: Yulong Kashi River; cold and arid area; reanalysis data; temporal and spatial distribution character- 
istics; SWAT model 


